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Abstract 

We present a systematic group space scan of discrete Abelian flavor symmetries for lepton 
mass models that produce nearly tribimaximal lepton mixing. In our models, small neutrino 
masses are generated by the type-I seesaw mechanism. The lepton mass matrices emerge 
from higher-dimension operators via the Froggatt-Nielsen mechanism and are predicted as 
powers of a single expansion parameter e that is of the order of the Cabibbo angle 6c — 0.2. 
We focus on solutions that can give close to tribimaximal lepton mixing with a very small 
reactor angle ^13 ^ and find several thousand explicit such models that provide an excellent 
fit to current neutrino data. The models are rather general in the sense that large leptonic 
mixings can come from the charged leptons and/or neutrinos. Moreover, in the neutrino 
sector, both left- and right-handed neutrinos can mix maximally. We also find a new relation 
6*13 < O(e^) for the reactor angle and a new sum rule 6*23 ~ f + e/"\/2 for the atmospheric 
angle, allowing the models to be tested in future neutrino oscillation experiments. 
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1 Introduction 



During the past decade, solar P, atmospheric [2], reactor [3], and accelerator [1] neutrino 
oscillation experiments, have very well established that neutrinos are massive. Since neutri- 
nos are massless in the standard model (SM), massive neutrinos signal physics beyond the 
SM. In fact, the smallness of neutrino masses ~ 10^^ . . . 10~^ eV can be naturally connected 
with grand unified theories (GUTs) [5], via the seesaw mechanism [3 E], in which the 
absolute neutrino mass scale becomes suppressed by an energy scale close to the GUT scale 
Mgut ~ 2 X IQi^GeV ^. 

Current neutrino oscillation data (for a recent global fit see Ref. [lU]) tells us that the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) lepton mixing matrix t/pMNS [ill] can be ap- 
proximated by the Harrison-Perkins-Scott (HPS) tribimaximal mixing matrix f/nps [12] as 



In f/npS) the solar angle 6*12 and the atmospheric angle ^23 are given by 612 = arctan(l/v^) 
and $23 = 7r/4, whereas the reactor angle ^13 vanishes, i.e., 613 = 0. The actually ob- 
served leptonic mixing angles in ?7pmns Kiay then be expressed in terms of deviations from 
tribimaximal mixing [131 [H] ^ "nearly" or "near" tribimaximal lepton mixing [15j. 

Many models have been proposed in the literature to reproduce tribimaximal leptonic 
mixing using non-Abelian flavor symmetries (for early models based on and examples 
using the double covering group of ^4, see Refs. [I6] and [E]). These models, however, have 
generally difficulties (for a discussion see Ref. [18]) to predict the observed fermion mass 
hierarchies and the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix Vckm [ISj 
(see Ref. [20] for models including quarks and Ref. [21] for unified models). An interesting 
connection between the quark and the lepton sector, on the other hand, is implied by the 
idea of quark-lepton complementarity (QLC) [22], which is the phenomenological observation 
that the measured solar mixing angle very accurately satisfies the relation 612 + Oq ~ vr/4, 
where Oq — 0.2 is the Cabibbo angle. QLC has been studied from many different points of 
view: as a correction to bimaximal mixing ^23j, together with sum rules [23], with stress on 
phenomenological aspects [25], in conjunction with parameterizations of t/pMNS in terms of 
9c [26] , with respect to statistical arguments [27] , by including renormalization group effects 
[28] , and in model building realizations p^. 

In Refs. [30l|3T], we have suggested a generalization of QLC to "extended QLC" (EQLC), 
where the mixing angles of the left- and right-handed leptons can assume any of the values 
^,e, e^,0. Here, e is of the order of the Cabibbo angle e ~ 0.2. By expressing also the 
lepton mass ratios as powers of e, we have derived in Ref. [31] for the CP-conserving case (a 
discussion of nonzero phases can be found in Ref. [32]) in total 1981 qualitatively distinct 
mass matrix textures for the charged leptons and neutrinos that lead to nearly tribimaximal 
neutrino mixing with a small reactor angle ^13 ~ 0. For these textures, the neutrino masses 
become small due to the canonical type-I seesaw mechanism (for a related approach see 
Ref. [33]). The matrix elements of these textures are in the flavor basis all expressed by 
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powers of e, which serves as a single small expansion parameter of the matrices. This 
suggests a model building interpretation of the textures in terms of flavor symmetries, e.g. 
via the Froggatt-Nielsen mechanism [31]. 

In this paper, we describe the systematic construction of several thousand explicit lep- 
ton mass models, in which nearly tribimaximal lepton mixing and the mass hierarchies of 
charged leptons and neutrinos emerge from products of discrete Abelian flavor symmetries. 
Our motivation is that models with Abelian flavor symmetries generally have the merit that 
they need only a very simple scalar sector to achieve the necessary flavor symmetry breaking. 
We perform a group space scan of products of discrete flavor symmetries using the results 
of EQLC from Ref. |31j and restrict to the case of real lepton mass matrices, i.e., we consider 
the CP-conserving case. In order to generate small neutrino masses, we assume only the 
canonical type-I seesaw mechanism. The hierarchical pattern of the lepton mass matrices 
results from higher- dimension operators that are produced by the Froggatt-Nielsen mecha- 
nism. Moreover, we are interested only in flavor symmetries that yield nearly tribimaximal 
lepton mixing with a very small reactor angle 6'i3 ~ 0. 

The paper is organized as follows: In Sec. we introduce the notation for the lepton 
masses and mixings. Next, in Sec. [3l we specify our discrete flavor symmetries and describe 
the generation of lepton mass terms via the Froggatt-Nielsen mechanism. Then, we outline 
in Sec. HJour approach to the group space scan of flavor symmetries. Our general results, a 
list of explicit example models, a new relation for the reactor angle, and our sum rules for 
the PMNS angles including a new sum rule for the atmospheric angle, are shown in Sec. O 
Finally, we present in Sec. [6] our summary and conclusions. 

2 Lepton Masses and Mixings 

We assume the SM with gauge group Gsm = 5'[/(3)c x SU{2)l x f/(l)y plus three right- 
handed neutrinos that generate small neutrino masses via the type-I seesaw mechanism [7]. 
The lepton Yukawa couplings and mass terms are 

Cy = -iYe),,H*i,e'j - {Yj,),yHl,v'. - h^MR),^v'rv''^ + h.c, (2) 

where li = (z/j, e,)^, ef, and uf, are the left-handed leptons, the right-handed charged lepton 
doublets, and the right-handed SM singlet neutrinos, and i = 1, 2, 3 is the generation index. 
Here, H is the SM Higgs doublet, and Yd are the Dirac Yukawa coupling matrices of the 
charged leptons and neutrinos, and Mji is the Majorana mass matrix of the right-handed 
neutrinos with entries of the order of the B — L breaking scale Mb^l ~ 10^^ GeV. After 
electroweak symmetry breaking, H develops a vacuum expectation value {H) ~ 10^ GeV, 
and the mass terms of the leptons become 

Cr^^s = -{Mehjeie^j - {MDhji^iU^ - ^{MR)iju^u] + h.c, (3) 

where Mi = {H)Y£ is the charged lepton and = {H)Y£) ~ 10^ GeV the Dirac neutrino 
mass matrix. After integrating out the right-handed neutrinos, the seesaw mechanism leads 
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to the effective Majorana neutrino mass matrix 



(4) 



with entries of the order 10 ^ eV in agreement with observation. The leptonic Dirac mass 
matrices and Mjj, and the Majorana mass matrices M^i and Meg are diagonahzed by 



D 



R 



M, 



eff 



diai 



where Ui, Ui/, Uc, Ud', Ur, and f/^, are unitary mixing matrices, whereas M^"^^, 



(5) 

diag 



and M^^^^, are diagonal mass matrices with positive entries. The mass eigenvalues of 



the charged leptons and neutrinos are given by Ml 



diag 



diag(m 



and Mfi"^ 



eff 



diag(mi, m2, ma), where mi,m2, and ma, are the first, second, and third neutrino mass 
eigenvalues. We can always write a mixing matrix Ux as a product of the form 



Ux Dx Ux Kx 1 



(6) 



where Ux is a CKM-like matrix that reads in the standard parameterization (we follow here 
throughout the conventions and definitions given in Ref. [5U] ) 



Ux 



( C12C13 

~'5l2C23 ~ Ci2S23Si3e" 
\ ■S12S23 — Ci2C23Si3e 



5^ 



C12C23 - 
~Ci2S23 



S12C13 
- ■5l2'523'5l3e 



i(5^ 



■Sl2C23Sl3e 



i5^ 



— i5^ 

■523C13 
C23C13 



(7) 



with Sij = sinOfj, Cij = cosO^j, and O^j G {6*^2)^13)^23} all in the first quadrant, i.e. 



6ij G [0,|], and 6^ G [0,27r]. In Eq. ([6]), and i^'j,. denote diagonal phase matrices that 
are Dx = diag(e"^i , e"^2 ^ e"^3) and Kx = diag(e^°i , e'"^ ^ 1), where the index x runs over 
X = i,i',D,D',R,u. The phases in Dx and Kx are all in the range [0,27r]. The PMNS 
matrix reads 

f^PMNS = UjUu = t/pMNS-^Maj) (8) 

where f/ is a CKM-like matrix parameterized as in Eq. ([71), and Kuaj = diag(e^'^\ e"^^, 1) 
contains the Majorana phases (pi and (j)2- The CKM-like matrix f/pMNS in Eq- dH]) is described 
by the solar angle 612, the reactor angle ^13, the atmospheric angle ^23) and the Dirac CP- 
phase 6, which we identify in the parameterization of Eq. ([71) as 9f^ — > On and 6^ — >• 6. 



3 Flavor Structure from Symmetries 

Let us next extend the SM gauge group to Gsm Gp, where Gp is a flavor symmetry. We 
assume that Gp is a. direct product of discrete Z„ symmetries, i.e., 

Gf = Zn,X Zn^X ■■■ X Znm, (9) 

where m is the number of Z„ factors and the rik {k = 1,2, ... ,m) may be different. We will 
denote by \Gf\ = H'^j^n/j the group order (i.e., the number of elements) of Gp. Under Gp, 
we assign to the leptons the charges 

e- - ip\,pl...,pin) =P\ ^i- {qlql---,qin) = - ir\,rl...,r'j =r\ (10) 
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where the jth entry in each row vector denotes the Z^. charge of the particle and i = 1, 2, 3 
is the generation index (see Sec. [2]). In the following, we choose a convention, where for each 
group Zn^. the charges are non-negative and lie in the range 

plqlrle {0,1,2,..., n,-l}. (11) 

The sum Xk + Vk of two charges is equivalent with + Uk mod nk and only determined 
modulo Hk- In order to spontaneously break the flavor symmetries, we assume for each factor 
Zn^, a single scalar fiavon field that carries a charge — 1 ~ ra^ — 1 under Z^^ but is a 
singlet under all other Zn^ with j ^ k. Moreover, the fn^ are Gsm singlets. 

When the fn,. acquire nonzero universal vacuum expectation values (fuk) — "^^ Tion- 
renormalizable lepton Yukawa couplings and mass terms of the form 

Cy = -(n-ie«^)(F/).,if*£.e,^-(n-ie^^)(F^),,ia2i7£,z/; 

-^(nr=ie^-)MB_L(n)..^r^l + h.c. (12) 

are generated by the Froggatt-Nielsen mechanism by integrating out heavy fermions with 
universal mass ~ Mp, where e = v/Mp ~ — 0.2 is of the order of the Cabibbo angle, and 

a!lj = min {p^ + g| mod Uk, -p\ - Qj mod Uk], 

= min {gf + r^^ modnfc, -gf - r| modnfc}, (13) 
c'^j = min {rf + modn^, — rf — r| mod Uk}, 

whereas F/, Y^, and Y^, are dimensionless order unity Yukawa couplings. The modulo 
function mod Uk in Eq. ( JT3l) is a consequence of the cyclic nature of the Zn symmetries, and 
the minimum takes into account that the higher dimension operators can be built from both 
and the complex conjugated fields f*^. The important point is that the Zn^ charges 
of the leptons determine a hierarchical pattern of the Yukawa coupling matrices and the 
right-handed Majorana neutrino mass matrix. 

We define a "texture" as the matrix collecting the leading order products of e for a 
certain Yukawa coupling or mass matrix, thereby ignoring the information on the order 
unity coefficients Y^, Y^, and Y^. The lepton textures are therefore the 3x3 matrices with 
matrix elements approximating {Mi)ij, {Mp,)ij, and {Mji)ij, as 

m,,^nt=,e<, (M,,),,^n^=/^., (M^,),,^n^=,64. (m) 

In what follows, we will, for a certain model call the set of the three textures defined in 
Eq. ([HD a "texture set" . 

Note that in our models, the Z„ fiavor symmetries are global, but it might be important 
to gauge them to survive quantum gravity corrections [35]. The cancellation of anomalies 
for our symmetries could, e.g. be achieved by considering suitable extra matter fields which 
is, however, beyond the scope of this paper (for a recent related discussion of anomalies and 
other phenomenology, see Ref. [36]). 
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4 Scanning Approach 



In this section, we will describe how we identify among the models introduced in Sec. [3] those 
which give nearly tribimaximal lepton mixing. First, we pick some flavor symmetry group 
Gp as in Eq. ([9]) and assign to the three generations of leptons e^,£j, and v^, all possible 
charge combinations under Gp according to Eqs. ffTOl) and ffTTl) . Then, we determine from 
the charge assignments the corresponding lepton textures for M^, M/j, and Mr, following 
Eq. (031) . Next, to flnd models that can give a good flt to nearly tribimaximal lepton mixing, 
we compare the textures found in Eq. f[T^ with the list of 1981 representative texture sets 
given in Ref. [21]. In Ref. ^SlJ, we have, based on assumptions of EQLC, determined flts 
of the order one Yukawa couplings Y/, Y^, and Y^, to reproduce nearly tribimaximal lepton 
mixing along with a charged lepton mass spectrum me : : m,- = : : 1 and a normal 
neutrino mass hierarchy nii : m2 : m^, = : e : 1 in perfect agreement with current neutrino 
data (at the Scr confldence level (CL)). 

In Ref. [H], the textures have been extracted in the basis where f/^/ = 1, i.e., where 
the rotations acting on the el are zero. Although these rotations do not show up in the 
observables, they are important for formulating our explicit models using flavor symmetries. 
We therefore include now in our considerations textures of the charged leptons that have 
been extracted in bases where Ui' can be nontrivial, i.e., where we can have Uii ^ 1, thereby 
leading to a multitude of new explicit representations for the charged lepton textures. In 
complete analogy with Ref. [31], we assume for the mixing angles and phases entering the 
matrices Uii all possible combinations that satisfy 

C<3,<3e{0,e^e,^}, (15) 

while the phases can take the values 5^\ipi^ip2 G {0,7r}, whereas we can always choose 
(/jg = a{ = ^2 = (for ^ deflnition of the notation, see Sec. [2]). This applies the hypothesis 
of EQLC to the such that the left- and right-handed lepton sectors are now all treated on 
the same footing. As a consequence, after employing the texture reduction from Refs. [301 131]. 
we flnally arrive at an enlarged "reference list" of 43278 qualitatively different lepton texture 
sets. We use this list as our reference to match onto flavor symmetry models for nearly 
tribimaximal lepton mixing. Any flavor charge assignment that yields textures contained in 
this reference list provides a valid model that allows an excellent flt to nearly tribimaximal 
lepton mixingj^ 

We can now impose extra assumptions on the properties of the textures in our reference 
list to search for interesting flavor symmetry models. For example, we will demand that 
none of the textures is completely anarchic (or democratic) and that the charged lepton 
textures contain (after factoring out common factors) at least one entry e" with n > 4 (to 
have sufficient structure in the texture). This reduces the above reference list to a reduced 
list of 17772 distinct texture sets. In the next section, we will use this reduced reference list 
to perform the group space scan of flavor symmetries. 

^Note that in our reference list (after factoring out a possible overall power of e) we set an entry e" equal 
to zero when n > 3 (for neutrinos only) or when n > 5 (for charged leptons) . This is different from Ref. |31j , 
where such an entry is set to zero when n > 3 for both neutrinos and charged leptons. For specific flavor 
models, however, we will always show the actual suppression factor e", as predicted by the flavor symmetry. 
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5 Results of the Group Space Scan 



Let us now present the results of the group space scan for flavor symmetries that produce 
nearly tribimaximal lepton mixing. We assume throughout the models introduced in Sec. [3] 
and choose as flavor symmetriecl Gp = I^fcLi^nfc for ^ = 1;2,3,4. A complete scan has 
been performed for groups up to order 40 (for m = 1), 45 (for m = 2), 30 (for m = 3), and 
24 (for m = 4), with < 9 for m > 1. Valid models are selected as described in Sec. Hlby 
matching the textures generated by the flavor charge assignments onto the reduced reference 
list of 17772 non-anarchic texture sets. In total, we find in the scan 6021 such models that 




20 40 60 80 100 200 400 600 800 1000 1200 1400 



n Models n Models 

Figure 1: Number of flavor models leading to nearly tribimaximal lepton mixing as a function 
of the flavor group Gp for increasing group order. In the left (right) panel we have 10 < 
\Gf\ < 24 (24 < \Gf\ < 45). 

reproduce 2093 texture sets. The distribution of the models for the different groups Gp is 
summarized in Fig. [T]If| All these models allow for an excellent fit to nearly tribimaximal 
neutrino mixing (at 3a CL and most of them actually at la CL ^31j) with a very small 

^For an application of related flavor groups see Ref. [37] , 

^Additionally, we have also included the number of valid models for Gf ^ Z2 Z4 x Z5. 
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reactor angle 6^13 < 1° and the lepton mass ratio^ 

rUe : : rUr = e'* : : 1, mi : 1712 : = : e : 1, (16) 

i.e., we have have a normal neutrino mass hierarchy with a ratio Amg/Am^^jj^ ~ of solar 
over atmospheric neutrino mass squared difference. Since the neutrinos have a normal mass 
hierarchy, renormalization group effects are negligible (for a discussion and references see 
Ref. [SU). 

In Fig. [H we can observe the trend that the number of valid models generally increases 
with the group order \Gf\ up to some periodic modulation. It is interesting to give a rough 
estimate on how large Gp has to be in order to more or less reproduce an arbitrary texture 
set. For this purpose, note that (after factoring out a possible common overall factor e") 
we restrict ourselves in Mp, and to the matrix entries {0,e^,e, 1}, and in to the 
entries {0, e^, e^, e^, e, 1}, i.e., we have 4^"'"^ • 6^ different possibilities for arbitrary texture 
sets. On the other hand, we have nine different charges per group, leading to \Gf\^ different 
possibihties for the charge assignments in (cf. Eq. f fTOj) ). Note that we use \Gf\ as a 
figure of merit: The larger \Gf\i the more possibilities for the charge assignments we have. 
In order to reproduce any texture set, we roughly estimate that the number of possibilities 
for the charges should exceed the number of possibilities for the texture sets, i.e., 

>4ir -6 -60. (17) 

This means that, for instance, four Zn factors with moderate n/c, such as Z2 x Z3 x Z4 x 
should be sufficient, or two Zn factors with large enough n^, such as Zi and higher. In Fig. [H 
we have much less possibilities on the l.h.s. where 10 < \Gf\ < 24, while we have on the 
r.h.s. 24 < \Gf\ < 45. In fact. Fig. [1] seems to suggest that Gf = Z^ x Zg with \Gf\ = 45 is 
already entering the regime estimated in Eq. (fT7|) . 

We have checked that for the 6021 valid models practically all (i.e., more than 99%) of 
the Yukawa coupling matrix elements of Y^',Y^, and Y^, (see Eq. (1121) ) lie in the interval 
between e and 1/e. With respect to the expansion parameter e of our models, these matrix 
elements can therefore indeed be viewed as order one coefficients. 

Out of our set of 6021 valid models, let us now consider a few examples. In Tab.[T], we show 
22 explicit valid models by listing the flavor group with complete flavor charge assignment 
and the resulting textures for M^, Mjj, and M/j. The rough guideline for choosing these 
22 models was to give one example for each texture set previously identified in EQLC pT] , 
with distinct Mr, a charged lepton mass spectrum as in Eq. (fT6!l . and "natural" order one 
Yukawa couplings Y^, Y'j-,, and Y^. The complete information on the corresponding mass and 
mixing parameters of the 22 models is summarized in Table [2] in the Appendix. 

Table [2] demonstrates that our models are very general in the sense that they can exhibit 
maximal mixings in the charged lepton and/or the neutrino sector, and that the left- and/or 
right-handed neutrinos can mix maximally. In addition, we can have maximal mixings in all 
sectors, not only between the 2nd and 3rd generation, but also between the 1st and the 2nd 
as well as between the 1st and the 3rd generation. 

''We are interested here in an SU{5) compatible fit. 
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All 22 models in Tables [T] and [2] lead to the following PMNS mixing angles in the range 

34° < ^^12 < 39°, 9rs<r, ^23 ~ 52°, (18) 



in agreement with neutrino oscillation data (at 3a CL). The models are thus characterized 
by a very small reactor angle close to zero, and a significant deviation of about +7° from 
maximal atmospheric mixing. 
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Z5 X Zg 


17 


e 


e3 1 


) 


•( 


IJ 


1 




/£2 £2 

£2 £ £2 


) 


(4, 3), (0,3), (4,0) 

(3,1), (2,0), (0,0) 
(2, 2), (2,1), (0,1) 


Z5 X Z7 


18 


( 


e3 e2 e2 

i; 






fe' e 

e 1 
e I) 


1 


e 


^2 £2^ 

£3 e £-^ 

,£2 £^ ly 




(4, 4), (1,2), (0,1) 

(1,0), (1,5), (0,6) 
(2, 3), (3,1), (0,3) 


Ze X Zj 


19 


( 


e2 e2 

i; 






c i ;) 




( 


£ £2 £'^\ 
£2 1 £M 

€> 1 1 




(0,1), (0,3), (4, 4) 
(2,1), (1,4), (1,2) 
(2, 3), (0,3), (0,0) 


Z5 X Zq 


20 


( 


£2 e2 e2 




•( 


e2 e 1 

e ly 




■( 


^2 e2 £2^ 

£2 £ £2 
,£2 £2 ly 


1 


(2.0) , (3, 5), (1,3) 
(0,4), (3,2), (4,3) 

(2.1) , (0,4), (2, 3) 


Z5 X Zq 
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21 


1 


e2 e2 £2 


1 










/e e2 e^ 
1 

e3 1, 


1 


(3, 4), (4, 4), (1,2) 
(3, 5), (4, 2), (4, 4) 
(2, 5), (1,3), (1,0) 


Zs X 


22 


1 


^2 ^2 ^3 

e ly 


1 


^1 


1 e 1 
Vl e3 1^ 


) 


■( 


^1 e3 i\ 

^3 ^ ^3 

,1 e3 i; 




(2, 6), (0,0), (0,1) 
(0,6), (1,1), (0,8) 
(1,0), (2, 5), (1,0) 





Table 1: A list of 22 valid flavor models for nearly tribimaximal lepton mixing. Shown are 
the explicit flavor charges under the flavor symmetry group Gp and the resulting textures. 
Possible common overall suppression factors have been factored out of the textures. 



Some of the models in Table [T] allow to extract very easily new sum rules for the lep- 
tonic mixing angles using an expansion in e. For example, we find for model 7^6 the sum rules 

while model #8 exhibits the relations 

e ,1 9, 9 ^ „ TT e ,1 5, 9 

For both the models #6 and #8, the solar angle satisfies the well-known QLC relation 
^12 ~ f ~ 0c/\/2. The reactor angle, on the other hand, becomes for these two models small 
due to an apparent suppression by a factor ~ 6q. Note that a similar suppression of the 
reactor angle has been found for the model in Ref. [38] • The atmospheric angle follows in 
both models #6 and #8 the new sum rule 6^23 ~ f + ^c/ v^j which predicts a deviation 
from maximal mixing by an amount of approximately Oq / \/2. This prediction makes these 
models testable in future neutrino oscillation experiments. For example, the deviation from 
maximal mixing can be established at 3(T CL by the T2K or NOi/A experiments [39j. In 
addition, one can measure the sign of the deviation from maximal mixing (the octant) with 
a neutrino factory at 3(T CL for sin^2^i3 > 10"^'^ or at 90% CL otherwise [10]. 

The relations for ^13 in Eqs. f|T9l) and fl20|) . however, are not stable under variations of 
the Yukawa couplings and must therefore be a result of exact cancellations between large 
contributions from different sectors. To render these relations for ^13 natural, it might, thus, 
be necessary to extend these models by non-Abelian discrete symmetries as demonstrated 
in Ref. [H]. 

Let us therefore have a look at model #5, which has the sum rules 

"-I-^I-P o^-^-o(^), ,21, 

The sum rules for 6*12 and 6*23 are to leading order just as in the two previous examples. The 
relation for 6'i3, however, is different from that in Eqs. (|T9l) and (J20i) : Now, 6'i3 is very small 
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due to a suppression by a factor ~ 6^. Unlike for models #6 and #8, the relation for 6*13 in 
Eq. fl2Tl) is stable under variations of the Yukawa couplings. This means that by varying the 
order one Yukawa couplings, Oi^ picks up only a small relative correction that is suppressed 
by a factor ~ 6q. The stability of the relation ^13 ^ C>{e^) under 1% variations of the order 



612 613 023 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 



Figure 2: Solar (left), reactor (middle), and atmospheric (right) mixing angles (in degrees), 
for model #5 in Tabled] as a function of the expansion parameter e. The points correspond 
to a random variation of the order one Yukawa couplings in F/, Y^, and Y^, by about 1%. 
The plot in the middle shows that the relation ^13 ~ 0{e^) in Eq. fl2T|) is not due to an 
accidental cancellation between large mixing angles. 

one Yukawa couplings is shown in Fig. [21 While 612 and 623 exhibit variations of the order 
(0.5 — 1)° in the limit e — 0, the variation of 613 remains <^ 0.1°. Nevertheless, the sum 
rules for 6*12 and 6^23 are in all three above examples in Eqs. (fT9l) . (!20|) . and (12T1) . in this sense 
stable, since the leading term in the respective expansions is it/ 4. 

In analogy with the groups, we have also performed a scan of ^7(1) symmetries. In 
order to see whether the cyclic character of the discrete groups is a relevant feature in the 
construction of valid models, we have compared single Z„ groups Gp = Zn-^ up to \Gf\ =40 
and product groups Gp = Zn^ x Zn^ up to IG^I = 18 with corresponding Abelian groups, i.e., 
with f/(l) (for Gf = ^„J and f/(l) xf/(l) (for Gp — Zn^ X Zn,^ by letting (the absolute value 
of) the respective individual U{X) charges vary in the whole range 0, 1, . . . , rifc — 1. Thereby, 
we have found 129 models for the cyclic groups and 24 models using ?7(1) symmetries. All 
of the 24 f/(l) models correspond to the group Gp = Z^ x Zq, for which we have obtained 
in total 21 models. However, even though the U{1) models produce textures that are very 
similar to those of model 7^6 in Table [1] (with varying and varying first row in Md), 
Mp, and Mr get in these examples highly suppressed by factors < e*. This analysis already 
indicates that the cyclic character of the flavor groups may be essential for facilitating the 
construction of realistic flavor models. 
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6 Summary and Conclusions 



In this paper, we have constructed several thousand exphcit flavor models for nearly tribimax- 
imal lepton mixing from products of Z„ flavor symmetries. In our models, small neutrino 
masses emerge only from the canonical type-I seesaw mechanism. Upon flavor symmetry 
breaking, the Froggatt-Nielsen mechanism produces hierarchical Yukawa coupling and mass 
matrix textures of the leptons. These textures are parameterized by powers of a single small 
symmetry breaking parameter e ~ — 0.2 that is of the order of the Cabibbo angle 9q and 
arises from integrating out heavy Froggatt-Nielsen messenger fermions. 

All flavor models realize the assumptions of EQLC and yield an excellent fit to nearly 
tribimaximal neutrino mixing with a very small reactor angle ^13 ~ 0. Moreover, they lead 
to the hierarchical charged lepton mass spectrum as well as to normal hierarchical neutrino 
masses. In our analysis, we have restricted ourselves to the most general CP-conserving case 
of real lepton mass matrices. We have performed a systematic scan of the group space for 
groups with up to four Z„ factors and a maximum group order of 45. As a consequence, we 
have found 6021 valid models that reproduce 2093 distinct texture sets. 

A characteristic property of our flavor models is that large leptonic mixings can come 
from the charged leptons and/or neutrinos. In the neutrino sector, maximal mixings can 
arise in both the Dirac mass matrix or in the heavy right-handed Majorana mass matrix. 
Generally, we can have maximal mixings between any two generations in any lepton sector. 

Among several explicit models, we have found a model that predicts a new relation 
6'i3 = (9(6'q), for the reactor angle. Moreover, we found models that all satisfy a new sum 
rule ^23 ~ f + for the atmospheric mixing angle, which makes these models testable 

in future neutrino oscillation experiments such as the T2K and NOz/A experiments or at a 
neutrino factory. 

We wish to point out that in this paper we have established a connection between a model 
building top-down approach using flavor symmetries and the phenomenological bottom-up 
approach of Ref . [31] . While Ref . [31] deals with the extraction of viable lepton mass textures 
that are in agreement with observation, the current work successfully matches the textures 
onto explicit flavor models, where the textures are predicted from flavor symmetries and 
their breaking. 

We believe that it would be interesting to study our sample of flavor models with respect 
to further model building aspects, anomaly cancellation, the inclusion of CP- violating phases, 
as well as in view of lepton flavor violation or leptogenesis. 

Appendix 

In this Appendix, we list the complete information on the mass and mixing parameters of 
the matrices that are generated by the models in Table [TJ Here, labels in both tables 
the same model. The data in Table [2] allows to fully reconstruct the exact form of the mass 
matrices of the 22 models following the notation of Sec. [2] (for further detailed examples on 
such reconstructions, see also Ref. [31j). 
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# 


mY /rnj:, 
mf/MB-L 




( a^' dt' Z)^' \ 
(^^12) C^13> C'23J 

((5^',Q;f ,Q;f) 


{ qD oD oD \ 
1^12 5 "l3' "23 J 

(5^,(^f,(^?,(^l>) 


f oD^ f\D' f\D' \ 
[^12 ) t>13 ) C^23 J 

(5^',<,0 


1^12' ^13' ^^23/ 

((5«, (^f , <ff, iff) 


1 


(e, 1, e) 


(6^6^6^) 

(7r,7r,7r) 


(f,e^o) 

(0,0,0) 


(o,6M) 

(0,0,0,0) 


(e, e, e^) 
(0,0,0) 


(tt, 0,7r, 0) 


2 




(0,0,0) 


(e,0,e2) 
(0,0,0) 


(0,0,0,0) 


(6,e^0) 
(0,7r,7r) 


(e,0,f) 
(0,0,0,0) 


3 




(e,0,e2) 
(0,0,0) 


(e,e,62) 

(0,0,0) 


(vr,0,0,7r) 


(e,e,f) 
(0,7r,7r) 


(e,e,f) 
(0,0,0,7r) 


4 


(e,e, 1) 
(e,l,l) 


(e,e,e2) 
(0,7r,7r) 


(i,e^o) 

(0,0,0) 


V 4 5 54/ 

(0,0, TT, 0) 


(0,e^f) 
(0,0,0) 


(0,e,e2) 
(0,0,7r,0) 


5 


(e, l,e) 
(e, 1, 1) 


(0,6,6^) 

(0,7r,0) 


(e,0,e) 
(0,0,0) 


(e,f,l) 
(7r,0,0,7r) 


(e^f,e^) 
(0,7r,7r) 


(0,f'0) 
(0,0, 0,7r) 


6 




(0, e, e^) 
(0,7r,0) 


(e, 0, e) 
(0,0,0) 


(7r,0,0,7r) 


(7r,0,0) 


(e, f ) 
(tt, 0,7r, 0) 


7 


(e^e,l) 
(e^e,l) 


(o,eM) 

(0,7r,7r) 


(e,0,e2) 
(0,0, tt) 


(0,0,0,7r) 


(0,0, tt) 


(7r,0,7r,7r) 


8 


(eM,e) 

(e^e,l) 


(0,7r,0) 


(e,e',0) 
(0,0,0) 


(vr,0,0,7r) 


(0,e^f) 
(0,0, vr) 


(7r,0,7r, 0) 


9 


(e, 

(e,l,l) 


(e,0,e2) 
(0,0,0) 


(e,e,0) 
(0,0,0) 


(vr,0,0,7r) 


(7r,0,7r) 


(i,e,0) 
(0,0,0,7r) 


10 


(e2,e,l) 
1) 


(6^6^62) 

(7r,7r,7r) 


(1,^2,0) 
(0,0,0) 


(- -) 
(0,0,0,7r) 


(6^0,62) 
(0,0,0) 


(e,e^O) 
(0,0,0,0) 


11 


(e, l,e) 


(7r,0,0) 


(f,0,e2) 
(0,0,0) 


(7r,0,0,7r) 


(e, e^, e) 
(7r,0,0) 


(e, e^, e) 
(tt, 0,0,0) 


12 




(0,e,e2) 
(0,7r,0) 


(e,0,0) 
(0,0,0) 


(e,f'i) 
(7r,0,0,7r) 


(e^f,e) 
(0,7r,0) 


(7r,0,7r,0) 


13 


(e^6,l) 
(e^e,l) 


(6^6^62) 

(7r,0,7r) 


(e,0,0) 
(0,0,0) 


V 4 I 14/ 

(vr,0,7r,0) 


(6,6^,0) 

(0,7r,0) 


(0,e^e2) 
(0,0,7r,7r) 


14 


(e^e,l) 


(7r,7r,7r) 


(e,0,e2) 
(0,0, tt) 


(f,e^o) 

(0,0,0,0) 


(e^I,0) 
(0,0, tt) 


(e^f,e) 
(0,0,7r,0) 


15 


(6^,6,1) 


(7r,7r,7r) 


(e,0,e2) 
(0,0,0) 


(1,0,6^) 

(0,0,0,7r) 


(e^6,62) 

(0,0,0) 


(e,e,e2) 
(tt, 0,0, tt) 


16 


(e, l,e) 


(6,6M) 
(7r,0,0) 


(e,0,e2) 
(0,0,0) 


(7r,0,0,7r) 


(f,o,o) 

(0,0, tt) 


(0,0,0,0) 


17 


(e^e,l) 
(e^e,l) 


(7r,7r,7r) 


(e,0,e2) 
(0,0,0) 


(1,0,6^) 
(0,0,0,7r) 


(e^o,62) 

(0,0,0) 


(e,0,e2) 
(0,0,0,7r) 


18 


(e^e,l) 

(e2,e,l) 


(6^e^62) 

(7r,0,7r) 


(e,0,e2) 
(0,0,0) 


V 4 ) 14/ 

(7r,0,7r,0) 


(e,e^e2) 
(tt, 0,7r) 


(0,e2,0) 
(0,0,0,0) 
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19 


/'c 1 1 ^ 


(7r,7r,7r) 


(1,0,6^) 

^"0 n n^i 
lU, U, Uj 


5454/ 
[IT, U, U, Uj 


(0,e,e) 

[0, U, TTj 


(e^6^f) 

[0, U, U, Uj 


20 




(0,vr,7r) 


(f,o,o) 

(0,0,0) 


145'= 54/ 
(vr,0,0,7r) 


(e^e^a 

(0,0,0) 


(vr,0,0,7r) 


21 




(7r,7r,7r) 


(f,0,0) 
(0,0,0) 


(0,6^1) 
(0,0,0,0) 


(0,0,0) 


(e^6,e2) 
(7r,0,7r, 0) 


22 


(6^6,l) 


(e^e^O) 
(0,0,0) 


(l,0,6) 

(0,0,0) 


(vr, 0,0,0) 


(0,7r,0) 


(e^f,0) 
(0,0,7r,0) 



Table 2: Supplementary information for the reconstruction of the mass matrices and Yukawa 
couplings of the models in Table [H Note that we have made use of the freedom to set 
= iff = ipD' = af = for i = 1, 2, 3 and j = 1, 2. 
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